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We study subexponential instability to characterize a dynamical instability of weak chaos. We show that a 
dynamical system with subexponential instability has an infinite invariant measure, and then we present the 
generalized Lyapunov exponent to characterize subexponential instability. 
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Chaotic dynamical systems with the zero Lya- 
punov exponent have attracted much interest in 
non-equilibrium statistical mechanics. Chaos in 
such dynamical systems is called weak chaos. We 
consider the Pomeau-Manneville map as an ex- 
ample of weak chaos, whose dynamical instability 
is not exponential one. Then, we define subexpo- 
nential instability using the average of the log- 
arithm of the separation of two nearby orbits. 
A dynamical instability of weak chaos is charac- 
terized as subexponential instability. We demon- 
strate the relation between subexponential insta- 
bility and infinite invariant measure. In particu- 
lar, we show that subexponential instability im- 
plies infinite invariant measure. The generalized 
Lyapunov exponent is studied to quantify a dy- 
namical instability of weak chaos. 



I. INTRODUCTION 

In the foundation of statistical mechanics, chaos plays 
an important role in deriving a stochastic description 
from microscopic dynamics^. The existence of an equi- 
librium state in dynamical systems requires ergodicity, 
i.e., time average being equal to space average. This is 
because macroscopic quantities, which result from time 
averages of microscopic observation functions, are almost 
equal to space average (constant). 

Non-equilibrium phenomena can be classified into two 
different states: non-equilibrium stationary state and 
non-equilibrium non-stationary state. In non-equilibrium 
stationary state, the time average of an observation func- 
tion is almost constant but it depends on subsystems. In 
other words, macroscopic observables restricted to sub- 
systems are almost constant but depend on the position 
of subsystems. On the other hand, macroscopic observ- 
ables are not constant in non-equilibrium non-stationary 
state. For example, an approach to an equilibrium state 
such as a diffusion is non-equilibrium non-stationary 
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state because macroscopic observables depend on time. 
Other examples are random processes such as turbulence 
and chemical reaction, where macroscopic observables do 
not approach an equilibrium state but change randomly 
in time. In such phenomena, the time average of an ob- 
servation function in underlying microscopic dynamics 
need to be intrinsically random. However, the founda- 
tion of non-equilibrium statistical mechanics on the basis 
of time average has not been studied at all. 

Recently, random behaviors of physical quantities 
in non-equilibrium non-stationary phenomena such as 
anomalous diffusions and intcrmittcncy have been stud- 
ied extensively in both experiments and theories^—. In 
particular, randomness of the generalized diffusion co- 
efficient has been clearly observed in the diffusion of 
mRNA molecules inside the cell^. Since then, this result 
has been elucidated theoretically through the continu- 
ous time random walk model of anomalous diffusion 4 . In 
general, anomalous diffusion processes can be represented 
by deterministic dynamical systems^. It has been known 
that dynamical systems describing subdiffusions, where 
mean square displacement increases as (x(t) 2 ) — o(t), 
show aging and infinite invariant measure 7 -. Infinite mea- 
sure systems, i.e., dynamical systems with infinite invari- 
ant measures, have attracted much interest to found the 
non-equilibrium statistical mechanics^ - — . The notable 
point in infinite measure systems is that the time aver- 
age of an observation function is intrinsically random^ii. 
More precisely, the time average of an observation func- 
tion converges in distribution, and its distribution de- 
pends on the class of the observation function. This is 
reminiscent of the randomness of the time average in non- 
equilibrium non-stationary states. 

In this paper, we develop the concept of ergodicity to- 
ward the foundation of non-equilibrium statistical me- 
chanics based on time averages in dynamical systems. 
Our aim is to characterize the dynamical instability of 
infinite measure systems. In finite measure systems a dy- 
namical instability can be characterized by the Lyapunov 
exponent. In a one-dimensional map T(x), the Lyapunov 
exponent defined by the time average of ln|T'(x)| indi- 
cates the exponential rate of the separation between two 
nearby orbits. By Birkhoff's individual ergodic theorem, 
the time average of ln|T'(x)| converges to the average 
of ln|T'(x)| with respect to the invariant measure for 
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almost all initial points^. Here, we focus on the subex- 
ponential instability which is characterized by the subex- 
ponential growth of the separation of two nearby orbits. 
Accordingly, the Lyapunov exponent converges to zero in 
dynamical systems with subexponential instability. We 
show that dynamical systems with subexponential insta- 
bility have infinite invariant measures. We then study 
the generalized Lyapunov exponent defined as the aver- 
age growth rate of subexponential separation. 



II. SUBEXPONENTIAL INSTABILITY 

Gaspard and Wang suggested that the separation be- 
tween nearby trajectories shows stretched exponential 
growth in Pomcau-Mannivelle maps with infinite invari- 
ant measures^. However, the dependence of the growth 
rate on an initial point was not studied. 

Based on the Darling and Kac theorem-^, Aaron- 
son proved that the time average of the positive L (to) 
function~ converges in distribution, where m is an invari- 
ant measure. That is, the growth rate strongly depends 
on an initial point and is intrinsically random. When 
the logarithm of the derivative of a map is the positive 
L^to) function, the normalized Lyapunov exponent con- 
verges in distribution: 



Pr ( — V In \T'(x k )\ <t\^ M a (t) 



asn-> oo, (1) 



where a n is called the return sequence determined by 
the map and M a (t) is the normalized Mittag-Lcfficr 
distribution^. Generally, the return sequence is deter- 
mined using the wandering rate defined by 



w n = m 



|J T- k A , 



(2) 



\k=0 



where < m{A) < oo. In particular, the return sequence 
is given by 

(3) 



r(l + a)r(2 - a)w n ' 

where w n is regularly varying at oo with index 

For example, the sequence a n of the Pomeau-Manneville 

map, 

X n+ l = T(x n ) = X n + CX l + p (mod 1), (4) 

where c > 0,p > 1, is given by 

np (p > 1) 



a n oc 



log n 



(P = l) 



(5) 



Because of ln|T'(a;)| € L l (rn), the distribution of the 
normalized Lyapunov exponent converges to a Mittag- 
Leffler distribution of order a, where the order a is given 




FIG. 1. Probability density function P a (t) of the normalized 
Lyapunov exponent of the Pomeau-Mannivelle map (p = 1.5). 
The sequence a n is given by a n = {\}n a , where a = 2/3 and 
(A) = 2.622. 



by I /p. Figure 1 shows that the distribution of the nor- 
malized Lyapunov exponent, A = — XX=o m l-^'C 2 "*)!) 
converges to a Mittag-Leffler distribution even when the 
time is finite (n = 10 8 ). 

We consider the separation of two nearby orbits {x n 
and x' n ) defined as 



A„ = lim 

Ai(0)->0 



Ax(n) 



Ax(0) 



n-l 



n r ^ 



k=0 



(0) 



where Ax(n) = x n — x' n . By |T]), the distribution for the 
growth of the separation can be given by 



Pr{A„ < t} -> M a 



bit 

a„. 



as n — > oo. 



(7) 



For large n, the average of the separation is given as 



(A n ) S (e a " 



I ML, 



(8) 



where (•) represents the average with respect to an initial 
ensemble, (-)ml represents the average with respect to a 
Mittag-Leffler distribution, and Y a is a random variable 
with a Mittag-Leffler distribution of order a. Accord- 
ingly, one can write 



(e 



We then have 



a n Y\ 



I ML 



y r(l + q) fc 

^ r(i + ka) ( 

fe=0 v ; 



(9) 



(e a - Y ) ML > J2 



k=0 



([ka\ +1)! 



(ry a )LpJ+i _ I 
> [ja^ 1 [( )bJ] (1Q) 

ja n - 1 

where 7 = T(l + a). The upper bound is obtained as 
1 n_1 

-lnA„ = -Vln|T'(x fc )|<M, (11) 
n n f — ' 

k=0 
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FIG. 2. (color online). Numerical simulations for the growth 
of the separation in Boole transformation (23). The sepa- 
rations \T'(x n ) ■ T'(x n -i) • ■ ■ T'(x\)\ for five different initial 
points are represented by five dotted lines, the logarithm of 
the average of A n and the average of the logarithm of A„ are 
represented by the thick line and the thick dotted line, and 
the theoretical curve Eq. (26) is also represented by the line. 



where M = ln[l + c(l + p)]. For large n, we have the 
following inequality: 



eMBn Ma ) < (A n ) < exp(Mn), 



(12) 



where B is constant. As shown in Fig. 2, the fluctuation 
of the separation A„ is extremely large and the average 
of A n is almost determined by the largest value of In A n . 
This implies that subexponential instability should be 
characterized as the average of the logarithm of the sepa- 
ration, (In A„), rather than that of the separation, (A n ). 
We call that a dynamical system has subexponential in- 
stability if there exists the sequence a n = o(n) such that 



/ 1 n-l > 
-5>|T'(z fc )l 



1 asn-> oo. 



(13) 



III. SUBEXPONENTIAL INSTABILITY IMPLIES 
INFINITE INVARIANT MEASURE 

We show that invariant measure can not be normal- 
ized when a dynamical system has subexponential in- 
stability, and then that the normalized Lyapunov expo- 
nent becomes intrinsically random. Here, we consider a 
one-dimensional conservative, ergodic, measure preserv- 
ing map with an invariant measure absolutely continuous 
with respect to Lebesgue measure in the sense of infinite 
ergodic theory^. In other words, we do not consider 
dynamical systems with the trivial zero Lyapunov ex- 
ponents such as T'(x) = 1 for almost all x. Moreover, 
we do not consider dynamical systems where there ex- 
ist attractive indifferent periodic points, (T m (x^))' = 1 
for period m points x*, ■ ■ ■ , x* m (k — 1, • • • , m) in the 



sense of Milnor 18 . The typical example of such dynam- 
ical systems is the dynamical system at the period dou- 
bling bifurcation points. At the critical point (m — > oo), 
(A n ) increases as power law-i^. However, we do not con- 
sider this case because there is no sensitive dependence 
on initial conditions nor an invariant measure absolutely 
continuous with respect to Lebesgue measure^. We note 
that there is at most one er-finite invariant measure in a 
one-dimensional conservative, ergodic, measure preserv- 
ing map. 

First, we assume Xk = \n.T'(xk) > 0. If X\, X 2l ■ ■ ■ are 
independent and identically distributed (i.i.d.) random 
variables, then by the laws of large numbers we have 



Pr 



■ x„ 



-E{X k 



> e> -> 



as n — > oo, 



(14) 

where E(X)(> 0) is the expectation of X, which con- 
tradicts subexponential instability (E(Xi + • • • + X n ) — 
o(n)). This implies that Xk is not an i.i.d. random vari- 
able. Moreover, by E(X\ + • • • + X n ) = o(n) we have 



E(X n ) 



as n — ¥ oo. 



(15) 



It follows that the probability Pi{Xk < e} increases with 
time for all e > 0. That is, x k is not an i.i.d. ran- 
dom variable, but trapped near the indifferent fixed point 
(T'(a) = 1) for a very long time^I, where a (< oo) is 
an indifferent fixed point. In other words, the empiri- 
cal distribution of an orbit {x k } converges to the delta 
distribution: For all d > b > a, 



p n {x)dx — > as n — > oo, 



where 



p n (x) 



S Xo (x) +S Xl (x) 



+ S Xn {x) 



(16) 



(17) 



Therefore, the invariant measure is an infinite one. This 
is because there exists a unique invariant measure abso- 
lutely continuous with respect to Lebesgue measure and 
the empirical distribution of an orbit converges to the in- 
variant measure if the map has a finite invariant measure 
(probability measure). 

As will be seen later, an indifferent fixed point of Boole 
transformation is at ±oo. In the case of T'(±oo) = 1, a 
dynamical system is an open system, i.e, the space of 
a dynamical system is not compact. By the condition 
of subexponential instability, E(X\ + • ■ ■ + X n ) = o(n), 
Eq. (16) is fulfilled for— oo<fr<<i<oo. Therefore, 
the invariant measure is an infinite one. This impossibil- 
ity of the normalization of the invariant measure stems 
from the non-compactness of the space of a dynamical 
system. Performing a change of a variable, we can ob- 
tain a compact dynamical system with indifferent fixed 
points^ (see Appendix A). Intuitively speaking, the in- 
finite invariant measure on the non-compact support is 
obtained by stretching the infinite invariant measure on 
the compact support near the indifferent fixed points. 
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Ni N 2 N 3 

FIG. 3. Schematic illustration of N n and M n 



Next, we consider the case of Xk E M., E(\Xi \ + ■ ■ ■ + 
\X n \) = 0{n) and E{X X + ■ ■ ■ + X n ) = o{rip-. Let 
Xi,X 2 , • • • be i.i.d. random variables with the mean 0. 
Then, there exists the sequence a n such that the limit 
distribution of (Xi + ■ ■ ■ + X n )/a n converges to be a 
normal or a stable distributional. It follows that the 
probability of the negative Lyapunov exponent (A < 0) 
remains positive and that the Lyapunov exponent con- 
verges to zero as n goes to infinity. This means that the 
dynamical system has an attractive indifferent periodic 
points. As we mentioned above, we do not consider this 
dynamical system because there is no sensitivity to initial 
conditions. Therefore, X\, X 2 , ■ ■ ■ are not i.i.d. random 
variables with the mean 0. 

There are two cases: (i) The expectation of X n is pos- 
itive and E{X n ) — > as n — > oo. In a similar way we can 
show that the invariant measure is infinite, (ii) Due to 
canceling out of Xx + • • • + X n , E(X\ + • • • + X n ) does 
not grow linearly. We introduce some notations (see Fig. 
3): the partial sum of X n , S n — X\ + ■ ■ ■ + X n , the 
maximum of S n , m n — max{0, Si, ■ ■ ■ , S n }, the num- 
ber of renewals of maximum, N n = minjfc : > 
mfc_i, mft > mk+i,k > N n -\}, and the nth maximum, 
M„ = m Nn . The condition E{X\ + • • • + X n ) = o(n) 
leads to E(N n - N n -i) = oo or E(M n - M n _i) as 
7i — > oo, where the interevent times, N n — iV n _i, are i.i.d. 
random variables. 

In the case of E(N n — N n -i) — oo and E(M n — 
M„_i) < oo, an orbit is occasionally trapped in the re- 
gion near the fixed point for a very long time because the 
residence time in the region near the fixed point becomes 
larger as the reinjection is closer to the fixed point. Since 
the probability density function of the residence time on 
[a — e, a + e] for the tail part is invariant for all e, the 
expectation of the residence time on [a — e, a + e] denoted 
by Ta € is infinite, where a is the fixed point. By the law 
of large numbers, we have 



n— 1 

i E 

n ^ 

k=0 



l[ a _ Cja+e ](x fc ) -> 



= 1, 



(18) 



for all e as n — » oo, where Ta c is the expectation of the 
residence time on the complement of the set [a — e, a + e] . 
This means that the invariant measure is an infinite one. 
It is noteworthy that the fixed point does not need to 
be the indifferent fixed point (see Fig. 4). As shown in 




FIG. 4. (Color online). The transformation with a flat critical 
point. The flat critical point is at x — 0.5 + 0. 



Appendix B, the transformation with a flat critical point 
such as 



T(x) 



2x see [0,1/2] 

ex P(-^T72+ 2 )/ 2 are (1/2,1]. 



(19) 



has subexponential instability with In |T'(x)| € K. and an 
infinite invariant measur o 25 ! 26 , which is compatible with 
the above result. 

In the case of E(M n — M„_i) — > as n — > 
oo and E{N n — N n -i) < oo, this leads to 
(T ( - N "- N ™- 1 ')(x Nri _ 1 )y — y 1 as n — y oo . It follows that 
an orbit is trapped near the indifferent periodic points. 
Hence, the invariant measure is an infinite one. Modify- 
ing the logistic map at the tangent bifurcation point to 
change the one-sided stable indifferent periodic points to 
the repulsive indifferent periodic points, we have a dy- 
namical system with an infinite invariant measure. In 
particular, the modified logistic map is given by 



T(x) 



-2(a - e) + (1 + 2\/2)e(l 
(1 + 2n/2>(1-x) 



e) x e [0, e] 



x e (e, 1], 

(20) 

where T 3 (e) = e (e < 0.2). The three iterated map T 3 (x) 
has three indifferent fixed points (Fig. 5). In particular, 



T 3 (x) 



C(x — e) + o((x — e) ) as x — > e + 0, where 



C > is constant. Therefore, the invariant measure is 
infinite one 2 ^. 

When In |T'(x)| is the positive L 1 (m) function, by the 
Darling-Kac and Aaronson theorem, the normalized Lya- 
punov exponent converges in distribution to a Mittag- 
Leffler distribution. On the other hand, when ln|T"(x)| 
is Lj oc (m) with finite mean function, the Lyapunov expo- 
nent converges in distribution to the generalized arcsine 
distribution 5 . When ln|X"(x)| is Lj (m) with infinite 
mean function, the Lyapunov exponent converges in dis- 
tribution to a stable distribution 5 -. However, these cases 
do not show subexponential instability. Therefore, the 
subexponential instability makes the Lyapunov exponent 
intrinsically random, and the normalized Lyapunov ex- 
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by 



FIG. 5. (Color online). The modified logistic map. The solid 
line and the dotted line represent T(x) and T 3 (x), respec- 
tively. 



ponent could possibly be described using a Mittag-Leffler 
distribution. 



IV. GENERALIZED LYAPUNOV EXPONENT 

Indicators characterizing weak chaos have been stud- 
ied by quantifying the complexity or the sensitivity to 
initial condition o 10 ' 27 i 28 . To characterize subexponen- 
tial instability we use the generalized Lyapunov exponent 
proposed by Korabel and Barkai 10 . The generalized Lya- 
punov exponent A(a„) is defined as 



X(a r , 



1 71-1 N 



(21) 



k=0 



where (•) represents the average with respect to initial 
points. It is remarkable that one can choose an arbi- 
trary initial ensemble if it is continuous with respect to 
Lebesgue measure. Note that there exists the sequence 
a n such that < A(a n ) < oo as n — > oo in dynamical 
systems with subexponential instability. The sequence 
a n presents a global instability of a dynamical instability 
and A = lim A(a n ) represents the rate of a subexponen- 

71— >00 

tial growth. We call (o„, A) the Lyapunov pair. 

We review the Darling-Kac and Aaronson theorem. 
Let T be a transformation on X with subexponential 
instability, and let ln|T"(x)| be the positive I/ 1 (m) func- 
tion. Then the normalized Lyapunov exponent converges 
in distribution to a Mittag-Leffler distribution: There ex- 
ists a n such that 

Prj — y^ln\T'(x k )\ <t\ ->M a (t) / \n\T'(x)\dm 

(22) 

asn^ oo. 

We demonstrate subexponential instability for Boole 
transformation. Boole transformation T : R — ¥ R is given 



X n J r \ — r(x n ) x 71 



(23) 



2n 



The invariant measure is uniform, i.e., of infinite mea- 
sure. The return sequence a n is given a&ii 



(24) 



By the Darling-Kac and Aaronson theorem, we can 
obtain the generalized Lyapunov exponent 

/oo 
In |1 + l/x 2 \dx = 2vr. (25) 
-oo 

Thus, the subexponential instability is represented by 

(In A„) ~ 2V2n as n -> oo. (26) 



V. CONCLUSIONS 

In conclusion, we have shown that dynamical systems 
with subexponential instability have infinite invariant 
measures. It follows that a dynamical system with a fi- 
nite invariant measure does not show subexponential in- 
stability. Considering Birkhoff's ergodic theorem, we see 
that a dynamical system has the positive Lyapunov expo- 
nent if the invariant measure m is a probability measure 
and the function In |T"(x)| is an L 1 (m) function. Further- 
more, we have presented that one dimensional maps with 
the non-trivial zero Lyapunov exponent are described as 
infinite measure systems such as transformations with in- 
different fixed or periodic points which are repellers and 
transformations with flat critical points. 

While it has been known that random dynamical sys- 
tems have aperiodic orbits with the negative Lyapunov 
exponent^, the relation between an invariant measure 
and the Lyapunov exponent is still unknown in random 
dynamical systems. However, the invariant measure be- 
comes infinite when the Lyapunov exponent is zero in 
the drift-diffusion model^S. Moreover, the transforma- 
tion Tx = Dx with probability 1/2 and Tx = x/D with 
probability 1/2 does not have attractors but have an infi- 
nite invariant measure absolutely continuous with respect 
to Lebesgue measure. We conjecture that the non-trivial 
zero Lyapunov exponent would imply an infinite invariant 
measure. 

Finally, we characterized the separation of nearby 
orbit as the generalized Lyapunov exponent. There 
arc different definitions of the generalized Lyapunov 
exponent^—. Kuzovkov et al. applied the general- 
ized Lyapunov exponent to the Anderson localization 
problem^. Compared to these definitions, the gener- 
alized Lyapunov exponent in this paper focuses on weak 
instabilities such as stretched exponential and algebraic 
instabilities. It would be interesting to make clear the 
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relation between the macroscopic observables in non- 
equilibrium state and the generalized Lyapunov expo- 
nent. We expect that the generalized Lyapunov expo- 
nent will be useful to understand the non-equilibrium 
non-stationary phenomena. In particular, we hope that 
random behaviors of transport coefficients in anomalous 
diffusions will be characterized by the generalized Lya- 
punov exponent. 



ACKNOWLEDGMENTS 

The authors wish to thank Tomoki Inoue, Takahito 
Mitsui, Masaki Nakagawa, Soya Shinkai and Nathan 
Weissburg for fruitful discussions. We would also like to 
thank Nickolay Korabel and Eli Barkai for reviewing the 
manuscript and discussing the results. This research was 
financially supported by the Waseda University Grant for 
Special Research Projects. 



Appendix A: Change of a variable 

In Boole transformation a change of a variable, 

<j)(x) = — — , yields a map with indifferent 

fixed points on [0, 1]. The transformed map S : [0, 1] —> 
[0,1], 



S(x) =(f>oTo<p- 1 (x), 



(Al) 



is given by 



x(l — x) 



S(x) = 



1 — x — x 

1 - 2x 
1 — 3x + x 2 



, x e [0, f/2] 



(A2) 



xe (1/2,1]. 



This map has two indifferent fixed points at x = and 
1, and the asymptotic behavior at x — is given by 



S{x) = x + x 3 + o(x 3 ). 



(A3) 



Hence, we see that the invariant measure is an infinite 
one. 



Appendix B: Subexponential instability of a map (19) 

Due to the flat critical point, the invariant measure is 
an infinite invariant measure and the invariant density is 
not bounded at x — 0. The inequality, 2~ JV ~ 1 < T(x) < 
2~ N , is fulfilled when the residence time in [0, 1/2] is N 
and x E (1/2, 1]. This inequality is solved as 



1 



iVln2 



1 

-<x< 



1 



(N- l)ln2 + 2 



(Bl) 



Using 



ln|T'(z)| 



In 2 



x G [0,1/2] 



L + 2-ln2 + ln(a;-i)- 2 x 6(1/2,1], 

X — n ^ 



we have the inequality for the partial sum Sn+i = X\ + 
• ■ ■ + Xn+i- 

21n{(iV- l)ln2 + 2} < S N +i = In \T'(x)\ + N\n2 
< -ln2 + 21n(iVln2 + 2). 

Let l n denote the number of reinjections on [0, 1/2] until 
n iterations. Since the distribution of the injection on 
[1/2, 1] is uniform, the distribution of the residence time 
in the interval [0, 1/2] is given as 



F(N) = Pv{N k <N} = Pr{T(.x) < 2^} 
Pr i x > 1 



Mn2 + 2 



1 - 



N\n2 



(B3) 



Therefore, 
E(N k ) - 



NF'(N)dN > 



1 



2 + ln2 



N 



: dN = oo. 
(B4) 

By the renewal theory^, there exists a constant K > 
such that 



In = K- + 

In n \mn 



as n — > oo. 



(B5) 



Using the arithmetic and geometric means inequality, we 
have 

E{X X + ■■■+ X n ) < E(S Nl+1 + ■■■ + S Nln+1 ) 

<2E[ln{(N 1 + l)---(Ni n + l)}] 
< 2E(l n In(n/Z„)) 



n ln ln n 

oc — as n 

ln n 



CO, 



(B6) 



where N x + 1 + • • • + N[ n + 1 = n. This means subexpo- 
nential instability, E(X\ + • • • + X n ) = o(n). 
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